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Abstract 

In this paper, we propose and analyze a spectrum sensing method based on cyclostationarity 
Oh I specifically targeted for receivers with multiple antennas. This detection method is used for determining 

q , the presence or absence of primary users in cognitive radio networks based on the eigenvalues of the 

cyclic covariance matrix of received signals. In particular, the cyclic correlation significance test is 
■ used to detect a specific signal-of-interest by exploiting knowledge of its cyclic frequencies. Analytical 

expressions for the probability of detection and probability of false-alarm under both spatially uncorre- 
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lated or spatially correlated noise are derived and verified by simulation. The detection performance in 
a Rayleigh flat-fading environment is found and verified through simulations. One of the advantages of 
the proposed method is that the detection threshold is shown to be independent of both the number of 
samples and the noise covariance, effectively eliminating the dependence on accurate noise estimation. 
The proposed method is also shown to provide higher detection probability and better robustness to noise 
uncertainty than existing multiple-antenna cyclostationary -based spectrum sensing algorithms under both 
AWGN as well as a quasi-static Rayleigh fading channel. 
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I. Introduction 

Spectrum sensing is a key step in effectively realizing cognitive radio networks (CRN). In the 
CR access paradigm, secondary users (SU) in a CRN are allowed to access spectrum reserved 
for use by licensed or primary users (PU) given that 1) those resources are either currently 
unoccupied or 2) interference to the primary network is kept under an acceptable level [2]. The 
main goal of spectrum sensing is to accurately and efficiently detect the presence or absence of 
a PU in a given band, usually under the constraint of a low signal-to-noise ratio (SNR). 

Several spectrum sensing methods have been proposed in the literature 0. In general, these 
methods can be categorized as being based on either energy detection, cyclic correlation (cy- 
clostationarity), or matched filtering. Energy detection requires the least prior knowledge about 
the signal, while matched filtering requires the most. Cyclic correlation-based techniques lie 
in between, requiring either prior knowledge or accurate estimation of the cyclic frequencies 
present in the PU transmission signal. Although energy detection offers the lowest computational 
complexity and is the optimal blind detector in the presence of i.i.d. noise, its performance relies 
on accurate knowledge of noise power due to the SNR wall phenomenon ||4). The detection 
performance of energy detection also degrades in a temporally correlated noise environment. 

In some scenarios, such as very low SNR regime or when signal selectivity is important, cyclic 
correlation-based methods offer several advantages over other spectrum sensing approaches. 
Unlike energy detection, they do not suffer from the SNR wall issue. These methods are also 
resilient to temporally correlated noise and enable signal- selective spectrum sensing where the 
presence of signals-of-interest (SOI) can be detected based on their unique cyclic features due 
to their modulation type, symbol rate, and carrier frequency 0. 

One issue encountered with all spectrum sensing methods is the effect of fading in the channel 
between the PU and SU. There is a decrease in the probability of detection whenever the channel 
is in a deep fade. This can be alleviated by exploiting spatial diversity either through the use 
of cooperative spectrum sensing [6] or, if available, the use of multiple antennas. As a result, 
spectrum sensing algorithms exploiting multiple antennas have received considerable interest 0, 

0. 

Algorithms that leverage the cyclostationarity property have been applied in the past for mul- 
tiple antenna receivers. In 0, the sum of the spectral correlation for each antenna was proposed. 
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Such methods are considered post-combining techniques since knowledge of the channel state 
information (CSI) is not exploited. On the other hand, pre-combining techniques which utilize an 
estimate of the CSI to varying degrees have been shown to have better performance in a random 
channel. A method based on equal gain combining (EGC) was investigated in IfTOll which uses 
phase offset estimates to align the raw samples from each antenna. The aligned signals are then 
summed before finding the spectral correlation. Finally, a blind maximal ratio combining (MRC) 
scheme was evaluated in 0TJ which utilized the singular value decomposition (SVD) to find an 
estimate of the CSI and applied MRC on the raw samples. 

In this paper we propose a spectrum sensing algorithm based on the cyclic correlation sig- 
nificance test (CCST) designed for use in a multiple antenna system which we refer to as 
Eigenvalue-Based Cyclo stationary Spectrum Sensing or EV-CSS. The CCST was used in lfT2ll to 
perform cyclostationary source enumeration using an information-theoretic criterion. However, 
the use of CCST in the context of multiple- antenna cyclostationary spectrum sensing has not 
been investigated in prior work. The performance of this method in fading channels has also not 
been evaluated nor compared to other spectrum sensing schemes that exploit cyclo stationarity. 
In this paper, we derive the analytical performance of this detection method in both AWGN and 
flat-fading channels. These expressions are then verified through simulations. The results also 
enable us to investigate sensing performance in a spatially correlated noise environment. 

The rest of the paper is organized as follows. The system model is introduced in Section HI 
including a brief discussion of cyclo stationarity. The proposed algorithm is detailed in Section 
Un] including analysis of its detection performance under both AWGN and flat fading. Numerical 
results for various scenarios are presented in Section [IV] Finally, the paper is concluded in 
Section [V] 

Notation: | A| and tr(A) denote the determinant and trace of square matrix A respectively. 
denotes the (i, j)th element of the matrix B and I is the identity matrix. The superscripts * and H 
denote the complex conjugate and the Hermitian (conjugate transpose) operations, respectively. 
Given two random vectors x and y, we define cov(x, y) = ^{xy^} — E{x}E{y H }. Given 
column vector x, diag{x} denotes a square matrix with elements of x along its main diagonal 
and zeros everywhere else. We will use the notation A/" c (m, S) to denote a proper (circularly 
symmetric) complex multivariate Gaussian distribution with mean m and co variance X. Finally 
we use the notation y = 0(g(x)) to indicate that there exists some finite real number b > 
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such that Hindoo \y/g{x)\ < b. 

II. Background and System Model 

A. Background on Cyclostationarity 

A signal is considered to be cyclostationary if its statistical properties are periodic. Equiva- 
lently, if the cyclic autocorrelation function, defined as: 

=i-L AT f*A* + d x " ( ( - d e ~' 2 ' a ' dt ' (1) 

J — T 

is non-zero with some r for at least one a ^ 0, the signal is said to exhibit second-order 
cyclostationary property with a referred to as the cyclic frequency. 

For example, in BPSK signals, cyclostationary features exist at a = ^- and at a = ±2f c + jr, 
where T b is the symbol period, f c is the carrier frequency, and k E Z. Detailed analysis of the 
cyclostationary features for various digital modulations can be found in 0. 

B. Signal Model and Assumptions 

We adopt a similar signal model as that used in ifTTI . The spectrum sensing problem is to 
decide between two hypotheses: Ho, where the signal is absent; and Hi, where it is present. 
The received signal samples under the two hypothesis are given respectively as follows: 

, , \v(n), Ho 
x(n) = I (2) 

\s(n)+r](n), Hi. 

The received signal, sampled at a rate of 1/T S , forms M streams coming from each antenna 
with N samples each. This received signal is defined as x(n) = [xi(n), X2(n), . . . , XM(n)] T . 
The received signal is the superposition of P signal sources (including both the SOI and any 

interferer) and can be expressed in vector form as 

p 

yj{n) =} j hj {n) ® Sj (n) + r) (n) , (3) 

i=i 

where cg> is the convolution operation over n and r/(n) is the receiver noise denoted by rf(n) = 
[r)i(n), r]2(n), . . . , r] M (n)] T , where every rji is a purely stationary Gaussian random process 
CR^(r) = for any a ^ 0) with variance of a^. For simplicity, we restrict that only one 
PU transmission, si(n), is considered a SOI and that it is cyclostationary with a unique cyclic 
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frequency a = a . The channel experienced by each of the P sources is given by hj(n) = 
[hji(n), hj 2 (n), . . . , hjM(n)] T , where hj k (n) is the channel between the jth source and the A;th 
antenna. We assume that the channel, although unknown to the receiver, stays constant over the 
spectrum sensing interval. Subsequently, we define the average signal-to-noise ratio to be 

A E{h H h\ 

SNR — j H { . (4) 
E [r] H r]} 

C. Spatially Correlated Noise Environments 

In the case of spatially correlated noise, which can happen when there is substantial ambient 
noise in the band, following |fT3l . we model rj(n) to have a co variance matrix given by R w = 
cov{t7, 77} where 

\ a ri, i — 3 

{Rr)T)}ij — E{r]fr}j\ — < (5) 

Thus with p s = 0, the covariance matrix simplifies to a v I giving spatially white noise, while 
p s = 1 gives fully correlated noise over all antennas. Varying degrees of partial correlation can 
be achieved by setting < p s < 1. 

III. Proposed Algorithm 

In this section, we describe the proposed method. We focus on a single cycle frequency 
detection, but this approach could be generalized to multi-cycle detection. 

A. Canonical Analysis 

The key idea of this detection algorithm is based on the theory of canonical analysis. As 
discussed in |fB|. and subsequently utilized in [fl2|. the number of common factors between two 
M x 1 time-series vectors x(n) and y(n) can be inferred from the rank of the matrix 

where we define R xy = cov(x, y) given the two random vectors x and y. 

Canonical analysis (see lH"5l for details) aims to find the relationships between two groups 
of variables in a data set. Given two random vectors x and y, of length m and n respectively, 
canonical analysis aims to find at most min{m,n} pairs of (uj,Vj) such that the correlation 
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between the linear combinations, Ui = ufx and V, = vfy, is maximized. An additional restric- 
tion is that Ui and Vi must be uncorrected with Uj and Vj for i ^ j. These linear combinations 
are referred to as canonical variates. The canonical variates are sorted in decreasing order of 
correlation such that the first canonical variates, U± and V\, have the highest correlation. The 
correlation coefficient, pi, between Ui and Vi is referred to as the ith canonical correlation. This 
procedure of finding and v, can be efficiently performed using a singular value decomposition 
(SVD) and the square of the canonical correlations can be found by finding the eigenvalues of 
©• 

In the context of cyclo stationary spectrum spectrum sensing, canonical analysis provides 
us with a very powerful tool to optimally combine the signals from M antennas and find 
the canonical correlations, resulting from up to M mutually uncorrelated cyclo stationary 
signals. This can be accomplished using the Cyclic Correlation Significance Test (CCST) |[T2| 
by performing canonical analysis on x(n) and x(n — r ) e -i 27rQnT » for a given lag r and cyclic 
frequency a. By finding the canonical correlations between these two sets of data, we are in 
effect measuring the maximum amount of cyclic correlation for all possible linear combinations 
of the signals coming from the M antennas. A threshold can then be applied on the combined 
Pi's to determine the presence or absence of the PU. Additionally, some cyclic frequencies, such 
as those located on a = ±2/ c for BPSK, only appear in the conjugate cyclic correlation. These 
can also be detected by instead performing the canonical analysis with x*(n — T ^ e -j 2 ^ anT > _ 

Prior to performing the detection, we pick the lag r that provides the best detection perfor- 
mance based on the modulation format used by the PU. This could be done off-line by performing 
the maximization, t = argmax r \R" a (r)\. 

B. Algorithm Description 

The steps of the algorithm are summarized as follows: 

1) Estimate the covariance matrix of size M x M 

JV-1-7B 

Rxx(r ) = — V x(n)x ff (n-r ). (7) 

N-l-T ^ 

2) Estimate the cyclic correlation matrix using a cyclic cross-correlogram at cyclic frequency 
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a and lag r , defined as 

jV-l-ro 

R- (r ) = $> (») x > " r °) e^™ ^- (8) 

n=0 

We will refer to the r -lag covariance matrices for both conventional and cyclic autocor- 
relation function simply as R xx and R xx from this point for the sake of brevity, since 
other r are not utilized by the proposed algorithm. The dependence on r will be indicated 
explicitly whenever necessary. The CCST is then calculated by finding the matrix 

R R XX R XX R XX R XX . (9) 

3) Find the eigenvalues, fi = [/xf , /x|, . . . , ji 2 M \ , of R. 

4) Compute the test statistic by combining the eigenvalues using 

M 



x±U(i-tf), (10) 



and finally calculate the test statistic: 



i=\ 



TJt = -m\n\. (11) 



The factor m = N — M — 1 is used to scale the test statistic so that its distribution is 
independent of the number of samples used [HH Sec. 8]. 
5) Decision: 7^ x 7, where 7 > is a threshold chosen to achieve constant false alarm 

rate (CFAR) which will be discussed in the following section. 
Note that all x H (n) can be replaced with x T (n) if the conjugate cyclic correlation matrix 
is needed. We refer to the version of the algorithm that uses x H (n) as the non-conjugate 
cyclic correlation significance test (NC-CCST) while the other is the conjugate cyclic correlation 
significance test (C-CCST). For the test statistic of each, we will use the notations 7^ x and 7^ x , 
respectively. 

C. Distribution Under Hq and Constant False-Alarm Rate 

Two key parameters are used to evaluate the performance of spectrum sensing algorithms. The 
detection probability or P D is the probability of being at Hi and accurately detecting the PU 
(Pd — P r (7^ x > 7 I "Hi))- On the other hand, the false alarm probability, Pfa, is the probability 
of being at H and mistakenly detecting a PU (Pfa — P r (7^ x > 7 | "Ho))- 
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It has been shown in |fT6l Sec. 8] that the limiting distribution (N — > oo) of the test statistic 
(ITTb for real and normally distributed random vectors approaches a \ 2 distribution with degree- 
of-freedom M 2 . Following a similar proof, it can also be shown that for zero mean, complex 
Gaussian random variables, the distribution is also \ 2 w i m degree-of-freedom M 2 when using 
the NC-CCST and M(M + 1) for the C-CCST. 

Based on the distribution of 7^ under Ho, the detection threshold 7 can be set to achieve a 
desired P F a by satisfying 

POD 

/ f x 2(x)dx = P FA , (12) 

where 

I M 2 NC-CCST 
k={ , (13) 

\M{M + 1) C-CCST 

and f x \{-) is the probability density function (pdf) of a x 2 random variable with degree-of- 
freedom k. 

These asymptotic distributions are verified to closely match simulation in Fig.CQfor N = 4000. 
Due to the scaling factor in (fTT|) . the distribution is independent of N. The empirical pdfs for 
two different a 2 values are also shown to demonstrate how the test statistic's distribution under 
Ho is independent of noise power. 

As introduced in Section Ill-Ci spatially correlated noise happens whenever is non- 
diagonal. This could also be interpreted as having a transformed noise vector 

rj = At] s.t. a ri AA R = R^^, (14) 

where A is an M x M matrix that determines the spatial correlation among antennas. In order 
to see the effect of correlated noise on the distribution of 7^ under Ho, it is helpful to use an 
alternate interpretation of canonical correlation as given in [17, Eqn. 2.1] such that 

R R a ° 

^XX XX 

R" a H T>H 



A 



(15) 



Rx 

using the covariance matrix estimates given in © and ([8]) and A as defined in (flOl) . Using (TT4l) . in 
conjunction with the multiplication property of determinants, and determinants of block diagonal 
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Fig. 1. Verification of the asymptotic distribution of the proposed test statistic (C-CCST) under Ho with different number of 
antennas (M = {2, 3, 4}). These plots show the accuracy of the analytical expression under N — 1000 number of samples per 
antenna (SNR=-10 dB). 



matrices, we find the value of A for correlated noise (A c ) to be 
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(17) 



Therefore the test statistic under Ho is invariant to any linear transformation on the noise 
measurements. As such the same expression for P F a as well as the threshold for maintaining 
constant P F a is applicable even for correlated noise environments. 



D. Distribution Under Hi and Probability of Detection 

In this section we derive the distribution of the proposed test statistic, 7^, under Hi. We 
begin by summarizing the prior work in statistics leading to the derivation of the complete non- 
null distribution of the canonical correlations. We show how this result is parameterized by the 
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canonical correlation, p, between the two signals being tested. In the context of cyclo stationary 
spectrum sensing, this corresponds to the value of the signal's cyclic autocorrelation, R" s = 
E [s(n)s*(n)e~ j27ra ° nTs ] . We then derive the canonical correlation, p, resulting from R^ s with a 
given a v , N, R w and h. Using this canonical correlation, the complete distribution of the test 
statistic under Hi is derived. Once the distribution is found, it allows us to find the theoretical 
probability of detection given as: 

P D = p r (T^ > 7 | Hi) = 1 - F TA \m (j) , (18) 

where Fj-a^ (•) is the cdf of the test statistic under "Hi which we find in the rest of this 
subsection. 

Several works in the past have contributed to deriving the non-null distribution of the test 
of significance of the canonical correlations originally proposed by Bartlett lfT8l . In 0~4]|, the 
approximate means and variances of the highest eigenvalue were found in the asymptotic case 
with only one non-zero eigenvalue (jii ^ 0, = 0,i > 1 ) and normality assumption. The 
distribution of the actual likelihood-ratio criteria, of the same form as (ITT]) , was found in iTTVll . 
However, both of these results break down in the case of local alternatives, which correspond 
to the cases when the distribution under "Hi is very close to the null hypothesis. In relation to 
the CCST, this corresponds to having a very-low SNR, which is clearly the case of interest in 
the spectrum sensing problem. 

Finally, the non-null distribution of the likelihood ratio criteria for covariance matrix under 
local alternatives were found in |[T9l . This criteria is used to test the independence between two 
multivariate random variables by combining the canonical correlation into a single test statistic. 
We use the same criteria to test for independence between the signal of interest and a frequency 
shifted version of itself. After the publication of this distribution, several works have focused 
on eliminating the normality assumption ll20l . EE]. However, in the case of spectrum sensing 
for CR, it is more likely to deal with detection under very-low SNR. If the noise is assumed 
AWGN and the a v > a s , the received signal, x(n), is approximately normal and the results for 
canonical analysis assuming normality can be used. 

Following 033 Thm. 4.1], the cdf of the the test statistic under local alternatives is found to 
be asymptotically distributed as a non-central chi-square. We have, 

F T ^\ Hl (x) = F x ,2(x,M(M + 1),5 2 ) + 0{m- 1 ), (19) 
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where F x v(x, d, 5 2 ) is the cdf of the non-central chi-square random variable with degree-of- 
freedom of d and non-centrality parameter 5 2 . The actual value of the non-centrality parameter 
can be found as 5 2 = tr(6 2 ), where = A/mdiag{pi, p 2 , . . . ,Pa/}- Where, p 2 for 1 < % < M 
are the eigenvalues of © using ensemble averages. Thus, they are solutions to the equation 



The true distribution deviates from the non-central chi-square with lower number of samples 
as indicated by the additional O^m^ 1 ) term in (TT9l) . The complete distribution up to the order 
C(m~ 2 ) can be found in lTT9l Thm. 4.1] which for the sake of brevity is no longer presented 
here. This more accurate distribution is an expansion based on non-central chi-square random 
variables of higher degrees-of-freedom with non-centrality parameters of tr(6 4 ) and tr(6 6 ). As 
such, a very accurate expression for F-j-a\ Hl ( x ) can be calculated with knowledge of 0. 

Let's assume, for simplicity, that only one signal of interest s(n) has cyclic frequency a. 
Therefore, pi = for i > 1. This assumption applies in almost all cases since two commu- 
nication signals will, with high probability, have different cyclic frequencies. However, it is 
straightforward to extend these results to the case of multiple signals with exactly the same a, 
since this simply corresponds to additional non-zero pi and the derivations presented here are 
still applicable. Thus, these scenarios can be treated theoretically as being single signal scenarios 
as long as the number of signals is less than or equal to the number of antennas M. 

The channel, h, is also assumed to be constant over one sensing period. Therefore, the 
expression for the distribution derived in this subsection is for a particular channel instance. 
Later we extend these results for flat-fading channels by integrating the distributions over the 
statistics of the fading channel. 

With these assumptions, we proceed with finding the values of pi as a function of the channel h, 
the signal of interest s{n) and the noise covariance R w The received signal can be expressed 
as x(n) = s(n)h + r){n). Using our initial assumption that the signal has unit power and is 
uncorrelated with the noise, the asymptotic zero-lag covariance matrix then becomes 




'XX 




(20) 



R 



•XX 



cov (x(n), x(n)) 



(21) 



hh H E [s(n)s*{n)] +R 



(22) 



(23) 
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On the other hand, recalling that noise has no cyclic features, the asymptotic cyclic cross 
covariance matrix can be found as 

= cov{x(ri),x(n)e- j27raonTs } (24) 
= hh H E [s{n)s*{n)e- j27raonT3 ] (25) 
= hh H iC (26) 

Subsequently the conjugate version, Rg S *, could also be used by replacing s*(n) with s(n). 

We again use the alternate interpretation in (fT5l) to find that the canonical correlation under 
Hi as 

ll A " lh H 2 W 

= 1! 11 1 ss \ (27) 

where the noise covariance matrix is R w = a rj AA H . The details of the derivation of (1271) are 
provided in Appendix [A] In the case of spatially uncorrected noise, we have A = I M so that 
the true correlation becomes 

P J^m, (28) 

which is only dependent on the 2-norm of the channel coefficient, the noise variance a v , and the 
cyclic correlation for the chosen a. A derivation of R" s is given in ll22l for various modulation 
schemes including BPSK, MSK, and QAM. 

A similar derivation can be done for C-CCST by replacing R^ s with R^ s , . Thus for a single 
SOI scenario we have 

e = V^[p,o,...,o]i M . (29) 

Which when substituted to (fT9l) gives us the complete distribution. In Fig. [2l we show that the 
theoretical results are in very close agreement to the simulation. 



E. Rayleigh Fading 

Now that we have the complete distribution for both the null and non-null hypothesis with 
AWGN parameterized by a particular instance of the channel h and the noise covariance R^, 
we can proceed to evaluate the performance of the proposed algorithm under a fiat-fading 
environment. 
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Fig. 2. Verification of the asymptotic distribution of the proposed test statistic (C-CCST) under Hi with different number of 
antennas (M = {2, 3, 4}). These plots show the accuracy of the analytical expression under N — 1000 number of samples per 
antenna (SNR=-10 dB). 

In particular we use a flat fading channel model where the channel vector, h, remains constant 
during the whole frame of samples used in detection. This is described using a channel vector for 
the ith frame as hj = [rie^ 1 ,r 2 e^ 2 , . . . ,r M e :j9M ] T , where r n is a Rayleigh distributed random 
variable of unit variance and 9 n is a uniformly distributed random variable in [0,27r]. 

We first assume spatially uncorrected noise (R vv = a v I). Based on (128b , the true correlation 
coefficient, which eventually determines the non-centrality parameter in (fT9b is a function of 
||h|| 2 = h^h. The distribution of h can be seen as a complex normal random vector of size 
M x 1 which results in ||h|| 2 ~ x\u- Combining ([19b and (128b . we find the new cdf of our test 
statistic under Rayleigh flat-fading and spatially uncorrelated noise as: 



becomes a generalized x 2 r - v - instead. 

Although a closed form expression is very difficult to derive for such an expression due to the 
presence of the non-central chi-square, there is still some insight to be gained by numerically 
integrating (130b to arrive at Pd and Pfa expressions under flat fading. 




(30) 



In the case of spatially correlated noise, the random variable fi — ||A 1 h|| 



h H (AA^ph 
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F. Comparison With Existing Approaches 

The algorithms for multiple-antenna spectrum sensing based on cyclostationarity that are 
currently in the literature can generally be classified into two categories. The simplest method to 
do this is to find the sum of the spectral correlation test statistic estimated individually from each 
antenna 0. We refer to this approach as SUM-MSDF (where MSDF means Modified Spectral 
Density Function). The MSDF is defined as the spectral correlation function (SCF) normalized 
by signal energy as discussed in ifTTI . 

Another existing approach is to sum the raw samples from each antenna and then perform a 
single spectral correlation test. However, we encounter a problem when the channel is not simply 
AWGN but instead has random fading. In this case, each antenna will have some unknown phase 
offset and attenuation. Thus, simply adding the raw samples non-coherently would decrease 
the probability of detection. This problem is addressed in ifTOll by first eliminating the phase 
rotation of signal samples coming from each antenna. An estimate of the relative phase difference 
between each antenna is calculated by finding both the cyclic correlation of one antenna chosen 
as reference (auto-spectral correlation function or auto-SCF) and the cross-cyclic correlation of 
every other antenna and the reference antenna. The phase difference can then be extracted from 
these two. We refer to this method in our comparisons as Equal Gain Combining (EGC). 

Finally, Maximal Ratio Combining (MRC) is used in [fTT|. Blind channel estimation is achieved 
by taking the vector corresponding to the highest singular value of ([8]) as an estimate of the 
channel, h. The raw samples from each antenna are combined using 

m = (3D 

||h|| 

The cyclic correlation test is then performed on the combined samples y(n). This method is 
called MSDF with blind maximal ratio combining or BMRC-MSDF. It was shown to outperform 
the other techniques but at the cost of additional complexity due to the channel estimation and 
combining. One issue with this approach is the fact that the cyclic correlation is calculated twice. 
The first is used to blindly estimate the channel and the second to perform the detection on the 
combined samples. In contrast, the method proposed in this paper only needs to perform the 
first part of BMRC-MSDF, finding the eigenvalues, and then uses the eigenvalues themselves to 
infer the presence or absence of the PU. 
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G. Advantages of the Proposed Algorithm 

As with other cyclostationarity-based spectrum sensing algorithms, one major advantage of the 
proposed method is its robustness to the noise uncertainty problem. Since the noise is assumed 
to be stationary and does not exhibit cyclostationarity at any a ^ 0, its cyclic correlation 
approaches zero as N — > oo. Thus, the effect of any error in the noise power estimate on the 
detection probability can be eliminated by taking more samples. However, in the interest of 
conserving power and arriving at a timely decision, both of which are high priority in the case 
of CR applications, we aim to minimize N needed to achieve a target This presents another, 
more subtle, issue related to noise uncertainty. 

In the non-asymptotic scenario, the methods based on the SCF (BMRC-MSDF, EGC and 
SUM-MSDF) under T-L have been shown to depend on both iV and the noise power a v ||22l|. 
Therefore, the proper detection threshold is still a function of the noise variance. By incorrectly 
specifying this threshold, the detector could be at the wrong point in the receiver operating 
characteristic (ROC) curve. Equivalently, the target CFAR cannot be achieved. However, as 
previously discussed and demonstrated in Fig. Q3 the proposed test statistic is independent of 
both cr^ and N. Consequently, the threshold 7 only needs to be chosen once for a given number 
of antennas M to guarantee CFAR. This property has been shown for other eigenvalue-based 
approaches Il23ll . It derives from the fact that noise power estimation is built-in to the algorithm. 

H. A Note on Complexity 

We provide an approximate complexity comparison of the proposed algorithm with the best 
performing existing algorithm (BMRC-MSDF) by taking number of complex multiplications 
required for each under the same number of samples N. Since the cyclic covariance operation 
and the SVD are common to both algorithms, they are not included in the analysis. 

Assuming the MSDF is calculated using an Appoint Fast Fourier Transform (FFT) it requires 
in the order of N\og 2 (Ns) multiplications. In addition, (M + 1)N multiplications are needed to 
perform the MRC and normalization. Finally, the correlation in frequency uses NN s /2 multipli- 
cations. Thus, the BMRC-MSDF approach performs in the order of N(log 2 (N s ) + N s /2 + M+l) 
multiplications without taking into account the SVD and the cyclic covariance. 

In comparison, the proposed EV-CSS method finds the conventional covariance in addition 
to an eigenvalue decomposition (EVD) and the same cyclic covariance as BMRC-MSDF, or in 
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the order of NM 2 multiplications. The operation R^R^x m ® ls essentially the solution to a 
generalized linear system which can be seen as an LU decomposition requiring approximately 
2M 3 /3 multiplications. Therefore, the EV-CSS approach requires in the order of NM 2 + 2M 3 /3 
multiplications in addition to the common operations with BMRC-MSDF. Since M is typically 
much less than both N and N s , there is overall a significant decrease in complexity with the 
proposed algorithm. For example, if we take N = 4000, Ns = 128, and M = 2, (same parameters 
used in QU), the BMRC-MSDF requires ~296K multiplications while EV-CSS needs only ~16K 
multiplications, without counting the common operations. 

IV. Numerical Results and Discussion 

In this section, simulation results are presented in order to compare the performance of the 
proposed algorithm with the various existing techniques discussed in Section IIII-FI In addition, 
theoretical plots are included to further verify the analytical expressions of the proposed method's 
performance. 

For these simulations we assume that only one PU has a feature at the chosen cyclic frequency, 
ceo- This PU is assumed to be transmitting a BPSK signal at a carrier frequency f c = 80 KHz 
with symbol period of 25 fxs. Each antenna of the SU is sampled at a rate f s = 320 kHz. For 
all algorithms, the same cyclic frequency located at a = 2f c is used. This cyclic feature is only 
present in the conjugate cyclic autocorrelation which means the C-CCST statistic is used. This 
feature is chosen because it is the highest magnitude feature among all cyclic frequencies. The 
maximum cyclic autocorrelation at this cyclic frequency is observed at r = which is the lag 
we will be using for all EV-CSS simulations. 

A. Threshold Selection for EV-CSS 

One key advantage of the proposed EV-CSS scheme over other spectrum sensing schemes 
is the simplicity of threshold selection to achieve CFAR. As discussed in Section IIII-CI the 
detection threshold, 7, is only dependent on M and not on N or a v . As a result of this, 
in a practical implementation of EV-CSS, the threshold can be pre-calculated using only the 
X 2 distribution parameterized by the number of antennas. As we will show in the following 
subsections, the theoretically determined threshold achieves the desired Pfa and matches very 
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well with simulations. Thus, only a single threshold needs to be stored. For the following 
simulations we set the CFAR to P FA = 0.1 unless otherwise stated. 

B. ROC and Detection Probability Versus SNR 

We first consider M = 2 antennas in the SU. The channel between the PU and each antenna 
of the SU, is modeled as a quasi-static Rayleigh fading channel with channel h as described in 
IIII-EI The fading is assumed to be frequency-flat and remains constant during the whole frame 
of N = 4000 samples per antenna used for detection. The noise in the antennas is assumed to 
be distributed as a zero-mean circularly symmetric complex-Gaussian, rj ~ jV c (0, <t v T). The PU 
signal energy is assumed to be unit energy and a rj is chosen to achieve an average SNR across 
antennas as defined in ©. 

Using these assumptions the ROC curves under SNR = — 10 dB and N = 4000 samples for the 
proposed algorithm and the other cyclic-based approaches are shown in Fig. [3] for comparison. 
Both theoretical results and Monte Carlo simulations (with 50,000 trials) are shown for EV-CSS, 
while only the Monte Carlo simulations are shown for the other methods since these are already 
analyzed in the respective works that proposed them. We have verified that these results agree 
with simulations presented in these prior work under similar assumptions. 

As seen in Fig. [3l there is very strong agreement in the theoretical and simulation results 
for the proposed algorithm. The performance of the proposed algorithm is clearly degraded 
when compared to a simple AWGN channel. However, it performs better than all the other 
techniques. Interestingly, the method also also outperforms BMRC-MSDF which, as discussed 
in Section IIII-Hl has significantly higher computational complexity. Although this result initially 
appears to be counter-intuitive, further experiments where only an AWGN channel is considered 
or if perfect channel state information (CSI) is assumed, show comparable performance between 
BMRC-MSDF and EV-CSS. Therefore, we can conclude that at very low SNR the blind channel 
estimates based on the SVD have large errors and the full benefit of MRC is not achieved. In 
contrast, the EV-CSS is able to fully take advantage of the information from all antennas because 
the algorithm works directly with the covariance matrices instead of utilizing an estimated CSI 
to pre-combine the signals. 

Another reason for the performance gain is that EV-CSS works directly with the time-domain 
cyclic autocorrelation of the PU signal which is particularly effective if some prior knowledge 
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Fig. 3. Receiver Operating Characteristic (ROC) of different cyclostationary-based spectrum sensing algorithms under Rayleigh 
flat-fading (SNR=-10 dB, N = 4000). 



about the cyclic frequency a is provided. An FFT based scheme requires the use of some form 
of frequency smoothing which degrades the cyclic feature. However, BMRC-MSDF can be made 
more robust to inaccurate knowledge of the cyclic features through these smoothing techniques 
as described in ifTTll . This has no bearing in the scenarios presented in this work since perfect 
knowledge of a is assumed and the parameters are chosen such that the cyclic frequency is 
perfectly aligned with an FFT bin. 

The EGC approach performs worst among all the techniques in a flat-fading environment 
because at very low SNR the estimation of the phase of h has very substantial error resulting in 
the test statistic being degraded in most cases. In contrast to this SUM-MSDF is able to separately 
calculate the spectral correlation of the signal on each antenna which offers significant gain after 
combining. 

The effect of varying SNR on probability of detection is also shown in Fig. [4} In both analysis 
and simulation plots for EV-CSS, 7 is set directly using (fT2l) in order to maintain a CFAR of 
P F a = 0.1. We again observe the strong agreement between theory and simulation for EV-CSS. 
In addition to this, the threshold selection is shown to be very effective in achieving the desired 
theoretical Pfa- As for the other techniques, the threshold needs to be determined as a function 
of SNR. This was achieved empirically in our simulations by Monte Carlo simulations of each 
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Fig. 4. Comparison of multiple antenna cyclostationary spectrum sensing techniques with varying SNR under Rayleigh flat- 
fading (M = 2, N = 1000, uncorrelated noise). 

algorithm under % . Although the distributions used for "Hi are designed for local alternatives 
(very low SNR), the analysis still matches simulation very accurately. This is true even at SNR 
> dB since the distribution converges to a non-zero mean Gaussian at very high non-centrality 
parameter which accurately describes the test statistic at both high N and high SNR. 

C. Varying Sample Size and Varying Number of Antennas 

The probability of detection over varying number of samples N E [1000, 5000] under an 
SNR of -10 dB and random flat-fading is shown in Fig. [5] In these simulations the Pfa = 0.1 
and the threshold is theoretically determined for EV-CSS. For the other techniques, the CFAR 
threshold is determined empirically through Monte Carlo simulatoins. As with varying SNR, 
only one threshold calculation is done for EV-CSS due to its independence to N. Significant 
improvements in detection probability of detection can be seen with increasing sample size. 
However, these gains tend to slow down with higher number of samples. 

The effect of number of antennas, M, on detection accuracy is studied in Fig. |6] In this 
figure only the best two algorithms (EV-CSS and BMRC-MSDF) are shown to facilitate the 
comparison. Note that for EV-CSS, to keep the the Pfa constant at 0.1 the threshold must be 
set to a new value based on (TT2)) . On the other hand, for BMRC-MSDF, the threshold is set for 
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Fig. 5. Comparison of multiple antenna cyclostationary spectrum sensing techniques with varying sample size TV under Rayleigh 
flat-fading channel (M = 2, N = 1000, SNR=-10 dB, uncorrected noise). 



different SNR and M. The a v across all antennas is assumed to be the same and an SNR is 
set using ©. Since both algorithms utilize some form of blind channel estimate we expect both 
to perform successively better as M is increased due to the improved spatial diversity provided 
by multiple antennas due to independent fading. Similar to previous results, the EV-CSS has 
better performance than BMRC-MSDF for different values of M at low SNRs. The agreement 
between theory and simulation for EV-CSS remains very strong even with higher M. 

D. Spatially Correlated Noise 

The effect of varying spatial correlation, p s , is shown in Fig. [7] In this simulation we again 
have M = 2 antennas. In the BMRC-MSDF simulations we assume perfect knowledge of the 
channel with h = [1, 1] T (no fading) and therefore the EVD-based blind channel estimation is no 
longer performed. This also results in BMRC-MSDF and SUM-MSDF having very comparable 
performance and thus, only one of them is shown. All methods, are degraded by increasing levels 
of spatial correlation. However, we see that BMRC-MSDF is more robust to such an impairment 
and ends up performing slightly better than EV-CSS at very high level of correlation. However, 
its more likely that each antenna would have only some slight level of correlation in a properly 
designed RF front-end, then the two methods could be regarded as being comparable in a simple 
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Fig. 6. The effect of the number of antennas on the detection probability of EV-CSS and BMRC-MSDF. The same number of 
samples per antenna N = 4000 is used. 
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Fig. 7. The effect a increasing spatial correlation on detection probability. A simple AWGN channel is used in these simulations 
to highlight the effect of spatial correlation (JV = 1000, SNR=-10 dB). 



AWGN channel. SUM-MSDF is seen to be the most robust to spatially correlated noise since 
it calculates the spectral correlation for each antenna individually. In fact, it shows some slight 
improvement at p s = 0.2. 
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E. Robustness to Noise Uncertainty 

Finally, we test the proposed algorithm's robustness under noise uncertainty and compare it 
with other cyclostationary-based spectrum sensing techniques. Noise uncertainty arises from 
the noise power level varying due to changes in the amount of thermal noise, amplifier gain, 
calibration error, and fluctuating interference [4J. The impact of the noise uncertainty is evaluated 
using the Bayesian statistics approach |]24]|. where a prior distribution on the noise power is 
assumed. In our simulations this is achieved by having an SNR with uniform distribution. 



where SNR is the average SNR which we set to be -10 dB. We then calculate the Pfa and 
P D by averaging over /snr(^)- The results for A G [0,3] dB are presented in Fig. [8] where 
the average Pfa = 0.1. As can be seen on the figure, both BMRC-MSDF and the proposed 
method perform at Pd = 1 when no noise uncertainty is present. However, BMRC-MSDF and 
EGC performance degrades significantly even with small A while EV-CSS is shown to be more 
robust to such an impairment. The slight degradation in EV-CSS is caused by averaging Pd over 
a uniformly distributed random variable SNR where the relationship between Pd and SNR is 
non-linear. As such, the values below SNR have a higher effect on the average Pd resulting in 
a net degradation. 



A multi-antenna cyclostationary-based spectrum sensing algorithm based on the cyclic corre- 
lation significance test was proposed and evaluated both analytically and through simulations. 
The method was shown to outperform existing multiple antenna signal- selective spectrum sensing 
methods in the literature. The computational complexity of the algorithm was also compared with 
that of the best performing existing algorithm that uses MRC by blindly estimating the CSI and 
was shown to require substantially less multiplications. The detection threshold for CFAR was 
also determined, both theoretically and via simulation, to be independent of the noise variance 
or the number of samples. This means that a single threshold is required for a given number of 
antenna, eliminating the need for separate noise estimation. The proposed method has also been 
shown to be highly robust to the effects of noise uncertainty. 




(32) 



V. Conclusion 
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Fig. 8. The effect of noise uncertainty, A, on the performance of various spectrum sensing schemes based on cyclostationarity. 
Noise uncertainty is assume to be uniformly distributed over an interval of [SNR — A, SNR + A] dB. SNR=-10 dB and 
N = 1000 samples (spatially uncorrected noise) 



Appendix 

Derivation of True Correlation under Hi 

In this appendix, we find the true correlation, p, under the non-null hypothesis Hi given a 
particular instance of the channel, h, and noise co variance matrix, R^^ = a v AA H . We begin by 
repeating (T23l) and (l26l) since these covariance matrices completely determine the test statistic. 

R xx = hh H + R vri (33) 

R«° = hh*iC (34) 

From (fl"5l) and under the assumption of only one SOI (pj = for i > 1) we have 
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where we recognize that both covariance matrices are Hermitian symmetric. Using the determi- 
nant for block 2x2 block matrices, we have 
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Substituting R xx and grouping together scalar terms 



Using Sylvester's determinant theorem for the sum of a full-rank and rank-1 matrix 
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To proceed further, we need to evaluate the inverse of (1331) which is a sum of a full-rank matrix 
and a rank-1 matrix. Using the results in ll25l Eqn. 1] the inverse becomes 
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Thus we have, 
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we arrive at the final expression 
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